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1. Introduction 



String theories present us with a natural opportunity to study extra dimensions. Within 
the Kaluza-Klein picture of the Universe these extra dimensions take the form of small, 
compact spaces with the full spacetime comprising four large dimensions with which 
we are acquainted and a compact manifold, Ai. While the large dimensions can be 
observed, the nature of the small manifold Ai cannot be directly seen, however some 
properties of M. can be deduced such as its dimension and, with certain assumptions, 
its topology. Provided there are no D-branes or fluxes, a situation we shall not con- 
sider in this paper, the manifold A4 must have special holonomy in order to preserve 
supersymmetry in the low energy theory|l|]. This restriction on the holonomy group 
implies that the manifold is a Calabi-Yau manifold. However there are a great many six 
dimensional Calabi-Yau manifolds and any one could be used for the compactification. 
Within each topology class there can be a further choice for the geometry, given by the 
value of the moduli parameters, where spaces with the same topology are grouped into 
moduli spaces and the continuous degeneracy of the moduh leads to massless fields in 
the low energy theory, moduli fields. Selecting one manifold to compactify on, given the 
huge range of possibilities, is a significant challenge for string theory; the observational 
consequences depend very heavily on JH, and without a unique manifold we lack the 
unique predictions we require. 

This situation was improved when it was realized that many of these Calabi-Yau 
manifolds were connected in moduli space 0, ^, ^, with the connecting geometry being 
singular. One way to picture these transitions between manifolds is to study the cycles 
within them, for example it may be that certain cycles collapse to zero size on one side of 
the transition and expand as different cycles on the other. Transitions may change the 
Hodge numbers, such as a conifold transition [§, || where an 5^ collapses and reappears 
as an or they may maintain the Hodge numbers but alter the linking numbers of 
the cycles, such as the flop transition [Q. Remarkably, string theory is able to make 
sense of these singular geometries by the appearance of new light states corresponding 
to D-branes wrapping the collapsing cycles [^. We can also investigate the low energy 
theory^] and deduce its dynamics]^, 0, 0, |T^, 0. The singular geometry has acted 
to connect the moduli spaces, its singularity is conical and takes the local description 
of a discrete quotient of a smooth manifold X, X = X/T, where F is a finite symmetry 
group; the singularity is then the fixed point set of F. Of particular interest are those 
singularities which allow a smooth resolution, by blowing up certain cycles of zero size 
contained within the fixed point set of F. This is achieved in practice by a surgery 
which replaces a ball around the conical singularity with a ball of a smooth special- 
holonomy space [|l^, |15|, [1^. The moduli of this special- holonomy space then appear 
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in the low-energy theory as moduh fields [O . 



We investigate the dynamics of the internal space while the cycles collapse, as in 
[pTSj we are interested in the gravitational effects surrounding the collapsing cycles, most 
notably whether a horizon forms to surround the cycle. The appearance of a horizon 
in the spacetime would indicate that the process acting to change the topology cannot 
be so easily understood as was previously thought, with the horizon acting to void 
the use of the currently-understood low energy theory close to the conical singularity. 
Previous studies have used the low energy dynamics of the moduli fields and inferred 
the results of the topology change from these; the formation of a horizon invalidates 
this approach. We will consider two geometries which act to change their topology 
by means of a conifold transition, namely the resolved and deformed conifods. We 
introduce these two spaces in the next section along with the nature of the transition 
we wish to instigate. We go on in sections ^ and ^ to describe the simulations we 
performed along with their results. We then attempt to ascertain the nature of the 
higher dimensional black hole in section H, and go on to give our conclusions in section 

i 



2. Conifolds. 

The conifold is a gravitational instanton exhibiting a singular point, the position of 
which we choose to be the origin. The utility of this space is that it represents the 



generic singularity within a class of Calabi-Yau gemetries |T9[ , and we know the explicit 
metric for the singular geometry [^] and its two smooth versions. These other two 
manifolds, called the deformed conifold and the resolved conifold, are similar to the 
singular conifold in that they have the same asymptotic structure, and they all admit 
a kahler metric. However, they differ from the conifold in that the singular point has 
been replaced in some way. The resolved conifold has an S"^ as a replacement for the 
tip, whereas the deformed conifold exhibits an 5*^. Following are the three metrics 
written within a specified coordinate system. 

2.1 Singular conifold 

A singular conifold is a Calabi-Yau threefold, it takes the form of a cone with a five- 
dimensional T^^'^^ base, 

ds^ = dr^ + r\dSbasef. (2.1) 
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Such a base has the form of S''^xS'^||^ and can be written in terms of two SU{2) 
groups, using two distinct sets of the conventional left-invariant one-forms of SU{2) 
which satisfy 

dcr* = -\eijk(r^ A CT^ 

dE* = -|e.,-fcS^' A T.K (2.2) 
With these we can write the base of the conifold 

d{s,asef = + al + + S^) + ha^ + Ss)^ (2.3) 

D 9 

2.2 Resolved conifold 

The resolved conifold is the alteration to the singular conifold achieved by expanding 
the singular point into an 5*^1211, It is further defined by a parameter determining 



the radius of the 5*^, we call this parameter a. In the limit that this radius drops to 
zero, we recover the singular conifold. Using the same one-forms to represent SU{2) 
the resolved conifold can be expressed as follows. 



+ l(6„^ + .=)p;-.E5)4(^)(.3-.E,)^. (2.4) 

At the origin this metric does not degenerate to a point, instead we find an S'^ 
which can be seen by setting r=0, 

ds'^\r=o = -dt^ + -dr^ + a\j:l + S^). (2.5) 
9 

This is a 2-sphere of area, Ana'^ (crossed with a 1+1 Minkowski geometry) which 
has replaced the origin of the singular conifold, and in doing so has smoothed the 
manifold. Clearly the resolved conifold approaches the singular conifold in either the 
high r or low a limits. 

2.2.1 Moduli space approximation 

In order to estimate the low energy dynamics of the resolved conifold we can allow the 
modulus a to vary in time to a small extent This involves introducing a new. 



time dependent modulus, which we called a{t) where a(0) = a. If we use the Einstein- 
Hilbert action and the resolved conifold metric, including this new dependence, we get 
an effective action: 

S.„=/A(|VS)^ (2.6) 
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This gives the approximation that a/ a{t) is hnear in time, implying that if we 
set a{t) on a course towards hitting zero, then the moduh approximation says that it 
will reach zero in finite time. As a = corresponds to the singular geometry then 
we see no obstruction to the geometry becoming singular within this approximation. 
By extending this approximation to include the full gravitational dynamics we hope to 
achieve a better understanding of this process. 

2.3 Deformed conifold 

There also exists the deformed conifold solution, which uses another method to rectify 
the singularity of the singular conifold |^, It involves expanding the singu- 

larity to the form of an 5*^. This is also described by a parameter, one which defines 
the radius of the S^, called e. 

Again using the conventional left-invariant one- forms which satisfy ( |!2.2D , and defin- 
ing K by 

^ ^ (sinh(2r)-2r)^ 7) 
23 sinh(r) 

we can write the deformed metric as: 



-df + 



3^3 



3K^ 



[0-3 + S3)' 



(2.8) 



+ ^sinh2 (D ((^1 + Si)' + (^2 - S^)^) 
+ ^cosh^ © ((^1 - Si)' + (^2 + S^)^) 



At the origin, this metric is also smooth as it does not degenerate to a point but 
to an 1^, 24|, this 3-sphere has replaced the origin and so removed the conical 
singularity. 

In the limit of high r values we can see that this tends to the same form as the 
singular conifold by noting that for large r 



(2.9) 



and so for large r. 



ds'^L- 



~dV + 



63 



2 



dr^ 



+ lei (2e-)^ Q {af + S? + + S^) + 1 (.3 + ^3)^) . (2.10) 
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This can be written, in terms of a new radial coordinate p, as 
.2 3 



P 



^et(2e-)^ (2.11) 



4 



ds\=^ = -dt" + dp" + p2 Q (^2 ^ ^2 _^ ^2 + S^) + 1 (as + Sa)'^ . 

This clearly shows that the asymptotic forms of both the singular conifold and the 
deformed conifold are the same, they only differ sufficiently close to the origin. 

2.4 Conifold transitions 

The deformed conifold and resolved conifold are each dependent on a continuous pa- 
rameter, which describes the size to which the new tip has been extended. The limiting 
case of these parameters gives the singular conifold. The common limit to these two 
manifolds allows for the transition (through the singular conifold) from a resolved coni- 
fold to the deformed conifold, despite their differing topologies p, 0; this is called the 
conifold transition. This transition can also occur in the opposite direction and it is 
these processes that we shall be investigating. 

We performed two classes of numerical simulations, one to simulate the collapse of 
a deformed conifold and one to simulate the collapse of the resolved conifold. Following 
from the results of these investigations we also construct a static black hole solution 
that we expect will be the final state of such collapses. 

2.5 Initial momentum 

For both cases of evolving a resolved or deformed conifold we know what the initial met- 
ric is, this leaves us with a choice for what the initial time derivative of the metric could 
be, its "momentum" . Noting that the spatial part of the metrics have co-dimension one 
- they depend only on r - we make the simplifying choice that the initial momentum 



also depends only on r. As pointed out in the fact that we are in higher dimensions 



means that such radial dependence still allows for gravitational waves, so the gravita- 
tional dynamics is not trivial. With this restriction of a radially dependent momentum 
the metric will remain co-dimension one. On top of this restriction we must also satisfy 
the constraint equations, known as the Hamiltonian and momentum constraints, given 
in the appendices. Our intention was to evolve the regular conifold geometries towards 
the singular conifold of section PTTI which we did by choosing a particular profile for 



the initial momentum, parametrized by an amplitude for the momentum, P, and a 
lengthscale over which the momentum was excited, Tq; the precise functional form is 
given in the relevant sections. These two parameters were the factors we varied in order 
to examine the collapse, in particluar we find the area of the black hole horizon as a 
function of these parameters. 
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3. The collapse of the resolved conifold 



To perform the simulations we need a suitable form for the metric which is general 
enough to give a consistent time evolution of the equations of motion, but with enough 
symmetry that numerical simulation is possible. For both the deformed and resolved 
conifold simulations we evolve a metric with four functions, each of which has an 
associated momentum. The equations of motion for these systems are given in the 



appendix and are written in such a way as to resemble the ADM formalism ||27|, p8 
3.1 Metric 

If we were able to solve the system analytically, a suitable form for the metric could be 
ds"^ = -dt^ + a^it, r)dr'^ + h^{t, r){al + al) 

+ c2(t,r)(S? + S^) (3.1) 
+ d\t,r){a^ + ^^f. 

However these may not be the optimum functions to evolve numerically. In par- 
ticular we note that if we take this form for the metric, then the initial conditions 
corresponding to the resolved conifold would mean that these profile functions con- 
tained singular regions, albeit co-ordinate singularities. In order to improve numerical 
stability, and to simplify the application of boundary conditions, we instead follow the 
method of [l^ and extract out various factors such that the functions we evolve are 
initially finite and asymptote to constant values, moreover they are symmetric under 

ds^ = -dt^ + A^{t,r)dr'^ + r'^B'^{t,r){al + al) 

+ {Qa^ + r^)C\t,r){T\ + T.l) (3.2) 

At r = we need to impose some boundary conditions, these follow by requiring 



local flatness at the origin |]26| and by maintaining A, C and D as even functions 
at the origin. 

A\t) = B\t) = C%t), (3.3) 
A{t,r) ~ A0(t) + O(r2), 

B{t,r) ^B\t) + 0{r^), (3.4) 
C(t,r) ~C°(t) + 0(r2), 
D{t,r) ~ D\t) + 0{r^). 
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3.2 Initial conditions 



By comparing ( p.2| ) and (|2.4| ), we can read off the initial conditions for the resolved 
conifold, 



A{0,r) = 

B{0,r) = ^, (3.5) 




C{0,r) 
D{0,r) 

This is a static metric, so if no initial momentum is added then no evolution occurs 
(this was used to test stability of our code). We added momentum and in doing so 
initiated the dynamical process, making sure that the momentum and Hamiltonian 
constraints were satisfied. Our specific algorithm was to add momentum to C (which 
describes the size of the two sphere) according to the form given below ( p.6| ). This 
depends upon two parameters giving the strength of the momentum (P) and the range 
the momentum extended away from the origin (ro). The momentum for B was taken as 
(|3.71) , in order to simplfy the initial conditions, this actually has the effect of maintaining 
a 3-sphere at the origin. 



C(0,r) = -P e (3.6) 

B{0,r) ^ DiO,r) 
5(0, r) D{0,r)' 

The form of the momentum is required to fall off sufficiently fast asymptotically 
so that we are not adding an infinite amount of energy, the exponential decay in ( p. 61 ) 
achieves that and means we could have a final state of a finite mass black hole. The other 
two initial conditions for the C and D momenta were determined by the momentum 
constraint and the Hamiltonian constraint, see the appendix for the equations. These 
initial conditions, along with the Einstein vacuum equations are enough to evolve the 
system forward in time. This was achieved in practice with a fourth order Runge-Kutta 
scheme, where the time step was changed dynamically such that: 

At = 0.005 lim A{t, r) Ar (3.8) 
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We also used a fixed spatial step of 0.004 with a total of 6400 grid points. Through- 
out the simulation we monitored both the Hamiltonian constraint and the momentum 
constraint, finding that they remained below 0.01. 

3.3 Apparent horizons 

We monitored the evolving system for the formation of an apparent horizon, a null 



trapped surface which implies the existence of an event horizon at some larger radius [p9|] . 
The apparent horizon, unlike the event horizon, can be detected locally and does not 
require knowledge of the full future evolution of the metric, it is defined as a congruence 
of null geodesies that no longer increase their area. For radial null geodesies we may 
write, 

d Area 







dt 
d Area 



dt 



null 



+ 



d Area 




dr 


dr 


t 


_dt_ 



(3.9) 
(3.10) 



null 



at the apparent horizon. This can be tested for at each individual timeslice and so 
reveal the existence of an apparent horizon and, if one is found, imply the presence of 
a horizon. We also measure the area of the apparent horizon at each time slice and so 
discover its evolution. For our simulations we find that its area increased monotonically 
but, as shown in Fig. |l], asymptotes to a constant value. This constant value we took 
to be a good approximation to the area of the resulting event horizon. 

3.4 Results 

We performed a number of simulations, varying the value of P and tq in the initial 
conditions of ( 3.6|) , examined whether an apparent horizon formed, and measured its 
area if one did; the results are presented in Fig. ^. For small values of P and tq we were 
not able to reliably extract a value for the final area of the apparent horizon and so this 
portion of the plot has been left blank. This was due to the large timescale needed to 
produce such black holes, and was beyond the dynamic range of our simulations. What 
we see from the data is that for all the cases which could be reliably tested we did 
observe the formation of a black hole. This is in contradistinction to the results of the 
system examined in |TB| for the evolution of the collapsing cycle in the Eguchi-Hanson 
geometry. 

In addition to the area of the horizon another property of the resultant black hole 
is the squashing of the angular part of the metric. This exists because the parameters 



relating to 6, c and d within (|3.1| ) need not be equal at the horizon, and their various 
ratios are referred to as the squashings of the metric at the horizon. The parameters 
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Time 



Figure 1: The area of the apparent horizon of the collapsing resolved conifold, changing in 
time. 




Figure 2: Resolved Conifold: The dependence of the final area of the horizon on the initial 
momentum profile. 

defining the extent of the squashing also change in time but, like the horizon area, 
converge as time goes on as shown in Fig. ^ for a particular example. The fact that 
both the apparent horizon area, and the squashing at the horizon are asymtoting to 
constant values is evidence that the final state is settling on a static black hole, we shall 
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present a candidate for these black holes in section ^. 



— (b/c)2 

— (c/d)2 

— {6/hf 



c 

(/) 



1.2- 



0.8 



°'^0 0.25 0.50 0.75 1.00 

Time 

Figure 3: The squashing values of (what was initially) the resolved conifold, seen to be 
converging in time. Note the slow change to the squashing (d/b) due to the initial condition 

(r 



3.4.1 Moduli space comparison 



In section \2.2.1\ we presented a prediction for the behaviour of a{t) based on the moduli 
space approximation, given that the static resolved conifold is a solution for any fixed 
a; this led us to conclude that ^ a{t) should evolve hnearly in time. We would like to 
check this result against the full numerical solutions that we have obtained, but this 
involves some ambiguity in defining what we mean by a once the metric has evolved 
away from the precise form of the resolved conifold. That is to say, once the profile 
functions have evolved away from the functional form given by ( |3.5| ), how do we extract 
a value for a? In practise we chose to define a{t) using the value of C(t, r) at the origin, 
C(t, 0) to extract a value for ^/a{t) by comparing expressions (|3.2| ) and ( p.4|) . 



y/^{t) = VQaoC{t,0). (3.11) 
As shown in Fig. ^ the time evolution of ^/a is indeed linear in the initial stages. 



4. The collapse of the deformed conifold 

In the previous sections we described the collapse of the resolved conifold, that was 
understood to corresdpond to the contraction of the 2-sphere which was made finite 
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Figure 4: The evolution of the the moduh of the resolved conifold in time. The approxima- 
tion predicts this shall be linear. 



in order to smooth out the conifold. Here we consider the second type of smoothing, 
namely the deformed conifold, which is achieved by having a finite size for the 3-sphere 
at the tip of the conifold. In the following sections we shall examine the time evolution 
of this geometry as this 3-sphere collapses. 

4.1 Metric 



In (p.8|) we presented the metric for the deformed conifold, so a natural choice for the 
metric that we should evolve is 



+ 6^ ((ai + Si)' + (as 
+ c' ((ai - Si)' + (^2 
+ (ag + Sg)'. 



(4.1) 



However, just as in the resolved case, this is not well-suited to a numerical solution 
owing to the profile functions diverging for the deformed conifold geometry. This is 
remedied by explicitly extracting the problematic terms, leaving us with smooth, finite 
functions to evolve, as well as allowing us to make the profiles symmetric under r —r. 
We therefore take the simulation metric to be 

ds"^ = -df + (^A{t,r) cosh (0)^^^^^ 
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+ (5(t,r)sinh(0) ((ai + Si)V(a2-S2)') (4.2) 
+ (C(t,r) cosh (0)' ((ai - Si)' + {a^ + S2)') 
+ (D(t,r) cosh (0)^(^3 + S3)' 

so that, y4, 5, C and D are all even at the origin and also they asymptote to constant 
values. This allows for more accurate application of boundary conditions arising due 
to maintaining even variables at the origin (implying (|4.3| )) and from requiring local 
flatness at the origin as described in ||26[ (leading to ( [4.4| )). 



A{t,r) ~ A°(t) + C(r2) 
B{t,r) ~ B\t) + 0{r^) 

D{t,r) ~ D°(t) + 0(r') 



(4.3) 



4.2 Initial conditions 

By comparing ( [4.2| ) and ( p.8| ), we can read off the initial conditions to be 



(4.4) 



A(0,r)2 
5(0, r)' 
C(0,r)' 
D(0,r)' 



4 



67^2 cosh' (I) 
Ij^ ^inh'(l) 
4 sinh' (I) 
1 4 cosh' (I) 

4 cosh' (0 

4 

6i^2 cosh' (I) ■ 



(4.5) 



As these are derived from the static deformed conifold these initial conditions will 
yield a static metric, as was indeed found in our simulations by way of a check for the 
code. We added momentum of the following form 



(7(0, r) _ Z)(0,r) 
C(0,r) ~ D{0,r) 



-P 



(4.6) 
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to initiate the dynamics, with the momenta for the A and B functions being determined 
by the constraint equations in the appendix. By choosing the C and D momenta to 
related we are making the choice of maintaining the form of the 3-sphere, at least 
initially. As in the resolved conifold initial conditions we note that the momentum 
decays away sufficiently fast that the energy input is finite, allowing for the possibility 
of forming a finite mass black hole. 



4.3 Results 



We used the same methods as described in section p.3| to search for the appearance 
of an apparent horizon, thus indicating the presence of a horizon. The results are 
qualitatively similar to the collapse of the resolved conifold, with an apparent horizon 
appearing in all the cases we examined. These apparent horizons had a time develop- 
ment similar to Fig. |^ and so we were able to establish a final value for their area in 
most cases, which we take to be a good approximation to area of the black hole hori- 
zon, this was checked by plotting the paths of outgoing null rays. The final area of the 
black hole depends upon the initial momentum we have added, which we characterized 
by a magnitude P and a spatial extent tq in ( |4.6|) ; varying these produced the graph 
displayed in Fig. ^ For small values of P or tq the area of the apparent horizon did 
not converge sufficiently rapidly to aquire an accurate value for the horizon area, so we 
have left such regions blank in Fig. ^ however, for all cases an apparent horizon was 
observed to form. 




Figure 5: Deformed Conifold: The effect on the final area of the horizon due to differing the 
initial momentum. 
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5. The resulting black hole 



For both the case of a collapsing resolved conifold and a collapsing deformed conifold we 
have seen that an apparent horizon always forms, for our choice of initial momentum, 
and this therefore implies the presence of a horizon. We have also seen that quantities 
such as the area of the apparent horizon, and the squashing modes of the angular part 
settle down to a constant value by the end of the simulation. This suggests that the 
final state of such collapses is a static black hole; in this section we suggest a possible 
candidate for the resultant black hole. 

We suggest a metric of the form 



-f{r)de + -—dr'' + h{r)\al + al) 



fir 



+ c(r)^(S? + S^) 
+ d{r)\a3 + S3) 



(5.1) 



as a candidate for the final state. So starting from this ansatz we found the equations 
for the profile functions by requiring Ricci flatness, these are given in the Appendix |y 
Unfortunately we were unable to derive an analytic solution to these equations, but 
we have been able to solve them numerically, showing that sensible solutions exist, one 
such solution with an event horizon at r = 1 is presented in Fig. At the location of 
the event horizon we have that / vanishes, with df/dr finite and non-zero, showing that 
it is simply a co-ordinate singularity. At the horizon we see that the profile functions b, 
c and d are all different, indicating that the horizon has a squashed angular geometry. 
We also see that the profile function, d, associated to the U(l) direction + S3 tends 
to a constant, rather than increasing linearly as b and c do. This shows that the spatial 
section of the black hole is not asymptotically a conifold. This behaviour was also 
seen in the case of collapsing cycles in the Eguchi-Hanson geometry ||18||, where a U(l) 
direction is picked out and aquires a constant radius asymptotically. 

These are candidates for the black holes formed by the collapse of the resolved 
conifold or, with appropriate coordinate transformations, the deformed conifold. We 
should of course be careful about the asymptotic structure of the solution which should 
not be altered by the collapse process near the origin. The discrepancy between the ini- 
tial condition being asymptotically conifold and the black hole asymptotics is resolved 
because the collapse takes an infinite amount of time to become the black hole, with 
a wave propogating outward from the collapse such that the interior is given by the 
black hole, and the exterior by the conifold. 
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Figure 6: The defining functions of the static black hole, with horizon conditions r = I, f 
4,/ = 0,6 = i\/3 and d= |. 

6. Conclusions 



Transitions capable of changing the topology of a Calabi-Yau manifold have been sug- 
gested previously [|], |], || and are a potential method for changing the topology of string 
vacua dynamically. These transitions would require the erstwhile distinct moduli spaces 
to share common, singular, manifolds. Dynamical application of the transitions would 
also require the cycles of the manifold to collapse successfully, that is to say without 
the formation of curvature singularities. Gravitational effects need to be accounted 
for as this collapse occurs, and failing to apply the effects of gravity may give false 
impressions for the transition process. In order to study the processes leading to the 
collapse of a two sphere or a three sphere, we have performed numerical simulations 
of two different, but related, cases of collapse. The resolved conifold, which involves 
a collapsing two sphere, and the deformed conifold, which involves a collapsing three 
sphere. These spaces are formed by two different methods of smoothing out a singular 
conifold, and it is using this singular conifold geometry that allows string theory to join 
the moduli spaces of these two regular, distinct geometries^. 

In this first study of the gravitational dynamics of conifolds, beyond the moduli 
space approximation, we have found that horizons are formed by the collapse of cycles 
in either of the two regular conifold solutions. These horizons take the form of two 
sphere crossed with a squashed three sphere. Such a horizon formed in all the cases we 



studied, which is different to the situation found in [IS] where a study of the collapsing 
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cycle in an Eguchi-Hanson geometry revealed cases where the cycle stopped coUpasing, 
with no horizon being generated. However, the case of the conifold allows for a much 
larger choice of initial conditions, and the ubiquitous formation of a horizon may be an 
artefact of our choice. 

Having established that horizons can form we have also found a possible candidate 
for the resulting black hole, by numerical means. This is a static black hole which also 
has a horizon comprising of a two sphere crossed with a squashed three sphere. An 
intriguing feature of this black hole, and one that it shares with the black associated to 
the collapse of the Eguchi-Hanson geometry, is that the asymtotic geometry contains 
a constant radius circle. 

In this paper we have only been interested in the gravitational dynamics coming 
from the Einstein-Hilbert action, and have not examined any stringy effects that may 
resolve the singularity. It is certainly true that if the scale at which the horizon forms 
is at the string scale, then these effects will be important. However, we have seen 
that a horizon forms for all of our initial conditions, so we should expect cases where 
the horizon size is greater than the string scale, and so field theory should be a valid 
approximation. For this initial study we have also neglected the dynamics from other 
supergravity fields such as the dilaton and form-fields. Indeed, this presents us with 
the natural extension of this work, namely to study the gravitational dynamics of the 
Klebanov-Strassler solution ||22|. This solution represents a supersymmetric deformed 
conifold that is held in place by the presence of fluxes, and has found many applications 
in the realm of flux compactification [^] and brane inflation |31|. Such calculations 



rely on the static Klebanov-Strassler geometry, without taking into account its own 
dynamics. If the dynamics of the Klebanov-Strassler geometry become singular then 
these calculations would need to be re-visited. 
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Appendices 



A. Resolved metric equations 

The components of the Ricci tensor for the resolved conifold case were calculated in 
the natural fiat basis, and act to define both the equations of motion and the constraints. 
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By defining the new functions 



b y 

Kb = T Db 



b " ab^ 
c d 

c cb 



d d' 
d db 

Ka = -, (A.l) 

a 

we find 

i?n = R22 = Kb + Kh {Ka + 2Kb + 2K, + K^) - ^ - (2D, + 2D, + D^) + 1 - 

D'. „ _ „ . 1 d^ 



i?33 = Ru = K, + K, {Ka + 2Kb + 2K, + Ka) - ^ - {2Db + 2D, + D. 



d 



a c 



2 



2c4' 

i?55 = Kd + Ka {Ka + 2Kb + 2K, + K^) - ^ - D^ (2D, + 2D, + D,) + ^ + ^, 

Err = Ka + Ka {Ka + 2Kb + 2K, + Krf) - 2^ - 2^ - ^ - (2^^ + 2D2 + Dj) , 

a a a 

Rtr = 2(^Db + KbDb^ + 2 (d, + K,D,^ + + K^D^ , 

Rtt = ~[Ka + Kl + 2kb + 2Kl + 2k, + 2Klka + Kf) . (A.2) 

Where the first four were used as evolution equations for the momenta, and the 
last two provided the constrint equations used to give consistent initial conditions, and 
to monitor the accuracy of the code. 



B. Deformed metric equations 

The components of the Ricci tensor for the deformed conifold case ( |4.1| ) were calcu- 
lated in the natural fiat basis, and act to define both the equations of motion and the 
constraints. By defining the new functions 

b b' 
Kb = - Db = -, 
b ab 

c r' 
K,= - D,= - 
c cb 

d db 

Ka = -, (B.l) 
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we find 

Ru = R22 = Kb + Kb {Ka + 2ir, + 2K, + K^, 
i?33 = ^44 = Kc + K, {Ka + 2K, + 2K, + ir^; 
i?55 = Kd + Kd {Ka + 27^5 + 2K, + iT^ 



Rtr 

R 



^-Dt {2D, + 2 
a 

^ - (2D5 + 2 
a 

^-Dd {2Db + 2 
a 



1 

262 + ; 


62 


C2 


Bc2d2 


862(^2 


1 

2c2 ^ ! 


C2 


62 


362^2 


8c2rf2 
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2d?^ 


d2 


62 


452^2 





-2 + - Ka)^ - 2 ( V + ^^^^^ " ~ ^ ~ 

= -{Ka + Kl + 2fi:, + 2Kl + 2E:, + 2i^2ierf + ifj) 

Where the first four were used as evolution equations for the momenta, and the 
last two provided the constrint equations used to give consistent initial conditions, and 
to monitor the accuracy of the code. 

C. Static black hole equations 

The need for Ricci fiatness imposes restriction upon the profile function of the metric 
(|5.1|) . These are calculated to be, 

/ 6' d d' 
f" = -V' T + - + W1 



b c 2d 



_ f'b' 4 d' b'b' 2b'c' b'd' 

~ T ^bf~ 2/63 ~ ~6 c d' 

f'c' 4 rf2 ^/^/ 2b' c' c'd' , 

^ =-V + -f-^ h T 

/ cj 2j& c b d 

f'd' d"" d'b' d3 
d" = -—^ + 2 — — - — + 



/ \Afb^ 6 4/c4 c 
_ '26' 2c' d!\ d^ d^ % % 
'T^T^7;^264 + 2^~62~^ 

6'\' fdV b'd b'd! d'd 
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